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[i]  A  global  ionospheric  electrostatic  potential  model,  which  we  refer  to  as  Nopper- 
Carovillano  (N-C),  can  be  linked  with  a  magnetospheric  potential  model.  The  latter 
model,  which  we  refer  to  as  Hill-Siscoe-Ober  (H-S-O),  computes  a  transpolar  potential 
$/>c(H-S-0)  based  on  solar  wind  parameters  and  region-1  field-aligned  currents  (FAC) 
from  the  magnetosheath  to  the  ionosphere.  The  polar  ionospheric  conductance  required 
by  H-S-0  is  defined  by  the  N-C  model.  In  this  way,  the  transpolar  potential  and  the 
associated  FAC  are  the  same  in  both  models.  A  distribution  of  region- 1  FAC  in  the  N-C 
model  predicts  a  two-cell  convection  pattern  which  is  in  reasonable  agreement  with 
plasma  drifts  measured  by  DMSP  (Defense  Meteorological  Satellite  Program)  satellites. 
The  H-S-0  model,  as  modified  by  N-C,  is  compared  with  the  Weimer  potential  model 
and  with  the  transpolar  potentials  observed  by  DMSP  satellites  during  the  6-7  April 
2000  magnetic  storm.  Good  agreement  is  found  in  both  cases.  The  region-2  (J2)  current 
is  estimated  from  the  Siscoe  (S-RC)  ring-current  circuit  model  which  is  driven  by  3>PC(H- 
S-O).  The  resistor  values  in  S-RC,  as  determined  by  N-C,  when  combined  with  the 
global  potential  solution,  make  it  possible  to  estimate  the  time  profile  of  the  equatorial 
penetration  electric  field  during  the  storm’s  main  phase.  With  the  values  obtained, 
shielding  occurs  within  1  hour  of  onset.  Equatorial  plasma  bubbles  (EPBs)  are  seen  some 
hours  after  the  initial  increase  of  $PC  and  are  qualitatively  consistent  with  the  equatorial 
penetration  electric  field  calculated  by  the  combined  model. 
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1.  Introduction 

[2]  Electrodynamic  coupling  between  the  interplanetary 
medium,  the  magnetosphere,  and  ionosphere  is  subtle  and 
complex.  In  this  paper  we  will  show  how  this  complexity 
can  be  approximated  by  joining  three  distinct  models  such 
that  penetration  electric  fields,  which  are  a  possible  source 
for  equatorial  plasma  bubbles,  may  be  estimated  from  the 
solar  wind.  Equatorial  plasma  bubbles  EPBs,  for  example, 
have  been  observed  in  both  ground-based  data  [Woodman, 
1970;  Fejer  et  at.,  1979;  Fejer  et  at.,  1991;  Whalen,  2001] 
and  satellite  measurements  [Fejer  et  at.,  1993;  Basil  et  at., 
2001 ;  Huang  et  at.,  2002],  Although  most  equatorial  plasma 
bubbles  are  generated  during  quiet  periods,  Greenspan  et  at. 
[1991]  and  Huang  et  at.  [2001]  noted  that  major  dismptions 
of  the  equatorial  ionosphere  also  occur  during  the  early 
phase  of  major  geomagnetic  storms.  EPBs  are  related  to  the 
low-latitude  ionospheric  electric  fields  which,  as  pointed  out 
by  Richmond  et  at.  [2003],  have  three  influences  of  com¬ 
parable  importance:  (1)  global  winds  driven  by  solar  heat¬ 
ing,  (2)  direct  penetration  of  polar  cap  electric  fields  to  the 
equator  that  are  partially  shielded  by  the  effects  of  region-2 
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(J2)  field-aligned  currents,  and  (3)  disturbance  winds  driven 
by  high-latitude  heating  and  ion -drag  acceleration.  In  this 
paper  we  focus  on  influence  2,  while  recognizing  the 
importance  of  influences  1  and  3. 

[3]  The  penetration  of  polar  cap  electric  fields  to  the 
ionospheric  equator  has  been  treated  by  Gonzales  et  al. 
[1979],  Spiro  et  al.  [1988],  Fejer  et  al.  [1990],  Fejer  and 
Scherliess  [1997],  and  Kelley  et  al.  [2003],  For  a  recent 
review,  see  Fejer  [2003],  Sophisticated  computer  simula¬ 
tions  by  Spiro  et  al.  [1988]  have  carefully  treated  the  time- 
dependent  magnetospheric-ionospheric  coupling  associated 
with  J2  according  to  the  insights  of  Vasyliunas  [1970]. 
Senior  and  Blanc  [1984]  developed  a  self-consistent  semi- 
analytic  approach  to  the  same  problem.  Scherliess  and 
Fejer  [1997]  used  the  Jicamarca  radar  observations  to  study 
the  characteristics  and  temporal  evolution  of  equatorial 
dynamo  zonal  electric  fields.  As  pointed  out  by  Richmond 
et  al.  [2003],  the  interplay  of  the  various  sources  for  the 
equatorial  zonal  electric  field  are  quite  complex  with  the 
ionospheric  winds,  at  times,  being  a  major  contributor  near 
dusk. 

[4]  The  purpose  of  the  present  report  is  to  provide  a 
quantitative  method  for  electrically  coupling  the  interplan¬ 
etary  medium  with  the  equatorial  ionosphere.  Our  method 
unites  three  methods  that  were  independently  developed. 
The  first  method  was  developed  by  Siscoe  et  al.  [2002] 
based  on  the  work  of  Hill  et  al.  [1976]  and  extended  by 
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Figure  1.  Data  for  the  magnetic  storm  on  6-7  April  2000.  (a)  The  dynamic  pressure  of  the  solar  wind 
plasma  at  the  LI  point,  (b)  The  interplanetary  electric  field  Esw  =  |Vx|  BT  sin2(y/2),  where  \\i  is  the 
“clock”  angle  between  By  and  Bz.  Figure  Id  shows  the  time  history  of  BT  =  \/(Bx2  +  By2), 
(c)  Finally,  the  ring  current  parameter  Dst  is  shown.  Also  shown  in  Figure  lc  are  the  times, 
denoted  by  open  triangles,  at  which  equatorial  plasma  depletion  bubbles  were  observed  during  the 
magnetic  storm.  The  solar  wind  variables  in  Figures  la,  lb,  and  Id  were  all  measured  at  the  LI 
point  and  appropriately  shifted  in  time. 


Ober  et  al.  [2003],  the  Hill-Siscoe-Ober  or  H-S-0  model. 
In  this  model,  the  potential  across  the  polar  ionosphere  $/>c 
and  the  associated  region- 1  field-aligned  currents  (FAC)  are 
determined  by  solar  wind  and  interplanetary  magnetic  field 
parameters.  The  second  approach,  developed  by  Nopper 
and  Carovillano  [1978]  (the  N-C  model),  detennines  the 
global  ionospheric  electric-potential  distributions  based  on 
the  region- 1  (Jl)  and  region-2  (J2)  currents.  The  third 
model  is  the  ring  current  circuit  model  (S-RC)  developed 
by  Siscoe  [1982],  In  this  model  (I>/>c  drives  J2  through  the 
ring  current,  which  is  presumed  to  have  inductive-like 
electrical  properties  consistent  with  the  observation  by 
Hines  [1963]  that  gradient  drifting  ions  are  gyroenergized 
by  the  cross-tail  electric  field.  We  join  these  three  models 
to  provide  a  time-dependent  method  of  estimating  the  level 
of  electric-field  penetration  of  the  solar  wind  to  the 
equatorial  ionosphere. 

[5]  The  penetration  of  interplanetary  electric  fields  to 
equatorial  latitudes  depends  on  the  relative  intensity  and 
location  of  J 1  and  J2.  In  short,  J 1  creates  an  eastward  electric 
field,  while  J2  creates  a  westward  (shielding)  electric  field  at 
equatorial  latitudes  near  the  dusk  terminator.  An  eastward 
electric  field  tends  to  enhance  the  Rayleigh-Taylor  instability 
associated  with  EPBs  (equatorial  plasma  bubbles)  while  J2 
tends  to  suppress  the  instability.  It  has  long  been  recognized 
that  region- 1  currents  are  closely  connected  with  the  trans¬ 
polar  potential,  while  region-2  currents  arise  from  the  injec¬ 
tion  of  plasma-sheet  ions  close  to  the  Earth.  The  effects  of 
high-latitude  changes  in  J 1  reach  equatorial  latitudes  at  MHD 
speeds  of  minutes.  On  the  other  hand,  the  timescale  for  the 
buildup  of  J2  is  on  the  order  of  hours. 


[6]  The  following  sections  briefly  review  the  H-S-0  and 
N-C  models  and  then  show  how  we  have  joined  them 
together  with  the  S-RC  model.  We  have  chosen  the  geo¬ 
magnetic  storm  of  6-7  April  2000  to  illustrate  the  com¬ 
bined  model.  This  storm  was  chosen  because,  occurring 
near  equinox,  it  allows  us  to  assume  symmetry  between  the 
northern  and  southern  ionospheres  which  therefore  removes 
complications  regarding  asymmetry.  We  demonstrate  that 
the  results  of  the  combined  model  are  in  general  agreement 
with  measurements  from  sensors  on  spacecraft  of  the 
Defense  Meteorological  Satellite  Program  (DMSP). 

2.  Observations 

[7]  Selected  geophysical  parameters  acquired  during  the 
6-7  April  2000  magnetic  storm  are  shown  in  Figure  1,  with 
the  solar  wind  parameters  shifted  in  time  from  the  LI  point 
by  the  x-component  of  the  solar  wind  velocity.  Figure  la 
shows  the  solar  wind  pressure  Psw,  Figure  lb  shows  the 
solar  wind  electric  field  Esw,  and  Figure  Id  shows  the 
component  BT  of  the  interplanetary  magnetic  field  that  is 
perpendicular  to  the  Earth’s  magnetic  dipole.  Dst  for  the 
same  period  is  shown  in  Figure  lc,  which  also  shows  the 
time  of  bubble  occurrence  during  the  stonn.  The  interplan¬ 
etary  electric  field  is  defined  by  Esw  =  |Vx|  BT  sin:’(\p/2), 
where  \p  is  the  “clock”  angle  between  By  and  Bz 
[Sonnerup,  1974]  and  Vx  is  the  x-component  in  GSM 
coordinates  of  the  solar  wind  velocity. 

[s]  Data  in  Figure  1  show  that  this  storm  period  was 
characterized  by  a  sudden,  strong,  and  sustained  increase  in 
Esw,  Psw,  and  therefore  in  $PC  (Esw,  Psw)  as  discussed 
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below.  The  sudden  increase  in  <J> PC  required  a  proportionate 
increase  in  the  portion  of  the  region- 1  FAC  directly 
connected  to  magnetospheric  boundary  layers  and  the 
magnetosheath.  The  region-2  system  grows  in  strength  over 
a  period  of  hours  due  to  a  delay  in  plasma  response  in  the 
inner  magnetosphere.  Note  that  the  electrical  equivalence  of 
inertia  is  inductance,  a  feature  used  in  S-RC.  The  storm 
followed  a  single,  sustained  pulse  of  energy  from  the 
interplanetary  medium  leading  to  an  immediate  and  simple 
main  phase  to  the  storm  as  indicated  by  the  Dst  trace. 

3.  Hill-Siscoe-Ober  Model 

[9]  Heppner  and  Maynard  [1987],  Weimer  [1995,  2001, 
2005],  and  others  have  described  the  distribution  of  poten¬ 
tial  across  the  high-latitude  ionosphere  based  on  measure¬ 
ments,  models,  theory,  or  some  combination  of  these.  One 
of  the  basic  components  of  any  of  these  is  the  transpolar 
potential  $PC,  the  difference  between  the  maximum  poten¬ 
tial  and  minimum  polar  potential.  Numerous  studies  exam¬ 
ined  the  response  of  (IVc  to  variations  of  the  interplanetary 
parameters  [e.g.,  Boyle  et  al.,  1997;  Ober  et  al.,  2003].  In 
this  paper  we  use  the  relationship  between  <&PC  and  the  solar 
wind  electric  field  Es^and  dynamic  pressure  described 
by  Hill  et  al.  [1976],  Siscoe  et  al.  [2002],  and  Ober  et  al. 
[2003],  Hill  et  al.  [1976]  originally  suggested  that  §PC  may 
be  represented  by  equation  (1). 


= 


<Fu<l)s 


(1) 


Here  is  the  transpolar  potential  that  would  result  from 
unhindered  dayside  merging  and  <J>5  is  the  saturation 
potential  obtained  when  region- 1  currents  distort  the 
dayside  magnetopause  to  inhibit  further  merging  [Siscoe 
et  al,  2002;  Ober  et  al.,  2003].  approaches  (bA/  during 
times  of  low  to  moderate  Esw.  <I> PC  approaches  $5  if  the 
Esw  becomes  very  large  such  as  during  a  major  geomag¬ 
netic  storm.  In  terms  of  the  solar  wind  parameters,  (1)  may 
be  rewritten  as  shown  in  (2)  using  definitions  described  by 
Siscoe  et  al.  [2002]  and  Ober  et  al  [2003],  Since  we  have 
included  the  viscous  interaction  term  proposed  by  Ober  et 
al.  [2003],  we  use  <l)PC(H-S-0)  to  denote  the  transpolar 
potential  (2) 


$PC(H-S-0) 


[30Pff  +  57.6E,w  ?£} 

[0.01 87  £,  £p  P£  +  0.036ESW  6,  £p  +  p£]  1 


[10]  The  units  of  Esw  and  Pw  are  mV/m  and  nPa, 
respectively,  and  Esw  has  been  previously  defined  above. 
In  (2),  £p  is  the  average  Pedersen  conductance  over  the 
polar  cap.  Ober  et  al.  [2003]  used  the  solar  radio  flux  at 
10.7  cm  (FI 0.7)  as  a  proxy  for  the  extreme  ultraviolet 
(EUV)  radiance  to  estimate  Zp  as  given  by  Robinson  and 
Vondrak  [1984],  In  the  present  work  we  use  the  Robinson- 
Vondrak  formula,  £p  =  0.88  (FI 0.7  cos  y)1/2,  with  the  solar 
zenith  angle  y  equal  to  80°  to  represent  the  average  polar 
cap  conductance  as  shown  in  (3). 


[11]  Given  a  set  of  values  for  Esw,  Psw,  and  F 1 0.7,  one  can 
use  (2)  and  (3)  to  estimate  the  magnitude  of  <FpC(H-S-0). 
Note  the  sensitivity  of  'IwfEI-S-O)  to  the  polar  cap  conduc¬ 
tance,  particularly  when  £p  is  on  the  order  of  a  few  Siemens. 
Siscoe  et  al  [2002]  use  Ohm’s  law  to  obtain  the  total  region- 
1  current  (Jl)  into  and  out  of  the  high-latitude  ionosphere  in 
response  to  the  applied  potential  drop  Pc(H-S-0).  For 
consistency  arguments  regarding  (2)  and  (4),  refer  to  Siscoe 
et  al.  [2002,  2004] 


Jl  =6,  E[  $pc(H  S  O) 


£  =  4.45  1.08  log  £[  (4) 


[12]  With  Siscoe  et  al.  [2002]  we  assume  (4)  is  the  total 
Jl  even  though  Jl  may  be  partially  observed  on  closed 
field  lines.  The  parameter  £,  scales  the  Pedersen  current 
closure  across  the  polar  cap  to  a  “single-wire”  approxi¬ 
mation.  Siscoe  et  al.  [2002]  did  not  try  to  describe  the 
distribution  of  current  into  and  out  of  the  region  1  system 
or  to  describe  the  resulting  distribution  of  ionospheric 
potentials. 


4.  Nopper-Carovillano  Model 

[13]  The  N-C  model  solves  the  current  continuity  equa¬ 
tion  to  determine  the  global  distributions  of  electric 
potential.  A  specification  of  the  distributions  of  Jl  and 
J2  FAC  and  global  ionospheric  conductance  is  required. 
Transpolar  potentials,  which  are  labeled  by  $PC(N~C)  in 
this  model,  are  determined  from  maxima  and  minima  of 
the  calculated  polar  potential  distributions  (see  Figure  2). 
Once  the  global  ionospheric  conductance  is  defined,  then 
the  N-C  model  determines  the  effective  polar  cap  conduc¬ 
tance  £,£p  by  finding  $PC(N-C)~]  for  Jl  =  1  MA,  J2  =  0, 
keeping  in  mind  that  $PC  is  in  kilovolts.  Consistency 
requires  that  this  value  for  £,£p  be  used  in  (2)  for 
determining  4Vc{H-S-0).  If  the  Jl  input  to  the  N-C  model 
is  set  equal  to  that  of  (4),  then  <I>Pc{N-C)  =  $/>c(H-S-0). 
The  joint  model  is  now  able  to  describe  the  penetration  of 
interplanetary  electric  fields  to  both  the  polar  and  low- 
latitude  ionosphere. 

[14]  The  N-C  model  provides  a  flexible  method  for  the 
distribution  of  the  total  field-aligned  current  (FAC)  flowing 
into  and  out  of  region  1  (Jl)  and  region  2  (J2).  The  specified 
Jl  and  J2  currents  enter  and  exit  the  ionosphere  but  are 
constrained  to  close  in  a  two-dimensional  conducting  shell. 
Siscoe  and  Maynard  [1991]  developed  a  similar  method  for 
distributing  FACs  and  deriving  ionospheric  potential  distri¬ 
butions. 

[15]  The  N-C  approach  uses  a  system  of  spherical  coor¬ 
dinates  whose  origin  is  at  the  center  of  the  Earth.  The  z-axis 
passes  through  the  magnetic  north  pole.  Initially,  a  three- 
dimensional  expression  for  J  is  reduced  to  two  dimensions 
by  imposing  the  condition  Jr  =  0.  Under  this  constraint  the 
two-dimensional  shell  current  J  and  conductance  tensor  £ij 
are  given  by 


£[  =  0.37i/F10.7 


(3)  Je  =  SeeEe  +  £&.pE..p  J9  =  — E&pEa  +  EWE,^  (5) 
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Polar  Cap  Potential  (kV) 


Figure  2.  A  comparison  of  the  Nopper-Carovillano  model  with  that  of  Weimer  [1995],  The  conditions 
are  for  equinox  using  the  default  N-C  conductance  model  (16).  The  N-C  results  for  the  polar  potential  are 
for  J 1  =  1.0  MA.  The  open  circles  show  the  solution  in  the  dusk-dawn  meridian  for  J2  equal  0.5  MA.  The 
open  triangles  are  for  zero  J2.  Note  that  the  effect  of  adding  J2  is  to  lower  the  polar  cap  potentials.  We 
have  included,  using  the  open  squares,  the  potentials  from  Figure  9  of  Weimer  [1995]  for  comparison. 
The  transpolar  potential  <&PC(N-C)  is  defined  to  be  the  be  the  difference  between  the  maximum  and 
minimum  voltages  and,  as  shown  in  the  text,  can  be  identified  with  the  transpolar  potential  of  Hill-Siscoe- 
Ober. 


where  Eo  and  E  .  represent  the  colatitude  and  azimuthal 
components  of  the  electric  field,  respectively.  The  con¬ 
ductance  related  tenns  are  defined  as 

EoEp  ^  EoEnsin(I) 

Lee  = . .  g - 

E„Epsin2(I)  +  (e2  +  E2  )  cos2  (I) 

= - V - - -  (6) 

D  =  Eo  sin2  (I)  +  Ep  cos2  (I) 

where  Ep  represents  the  Pedersen  and  £H  =  2Ep  represents 
the  Hall  conductance.  E„  denotes  the  conductance  parallel 
to  B,  which  Nopper  [1978]  takes  be  31.6  times  Ep.  The 
symbol  I  denotes  the  dip  angle  of  the  Earth’s  magnetic 
dipole  field.  We  presently  assume  that  the  dipole  axis  is 
aligned  with  the  Earth’s  spin  axis  and  perpendicular  to  the 
Earth-Sun  line.  Otherwise,  the  dip  angle  is  a  function  of 
geographic  longitude,  which  significantly  increases  the 
complexity  of  (6).  Therefore  as  the  model  now  stands, 
longitudinal  effects  are  not  included.  With  the  substitution 
E  =  — V<F  the  current  continuity  equation  reduces  to  a 
Poisson-type  equation  [Vasyliunas,  1970] 

V  •  [E  •  V4>]  =  j|  sin(I),  (7) 

where  j  ||  is  the  specified  distribution  of  FAC  associated  with 
J1  and  J2.  Note  that  jy  is  equivalent  to  a  positive  charge 


density  for  downward  current  and  a  negative  charge  density 
for  upward  current. 

[16]  Equation  (7)  is  solved  by  numerical  computation 
using  the  multigrid  method  [Press  et  al.,  1992;  Briggs  et 
al.,  2000;  Trottenberg  et  al.,  2001],  a  method  that  capitalizes 
on  the  fact  that  the  Fourier  components  of  the  potential 
converge  at  different  rates.  A  129  x  129  (16,641  points) 
global  grid  was  used  to  solve  for  the  distribution  of  iono¬ 
spheric  potentials.  Following  the  example  given  by  Nopper 
and  Carovillano  [1978,  1979],  the  J1  and  J2  currents  are 
distributed  in  two  crescent  regions  extending  120°  in 
magnetic  local  time  at  72°  and  66°  magnetic  latitude, 
respectively.  Each  crescent  region  is  centered  on  the  dusk- 
dawn  meridian.  Each  current  region  is  3°  wide.  For  our 
computations,  the  current  region  dimensions  are  32  by  2 
grid  spacings.  The  total  current  for  each  current  region  is 
uniformly  distributed  throughout  the  region.  In  this  way  a 
baseline  model  may  be  constructed  for  comparison  with 
data. 

[17]  Nopper  and  Carovillano  [1978]  used  a  globally 
varying  conductance  pattern  which  varied  from  a  solar- 
zenith  angle  of  -90°  to  +90°  with  enhancements  at  the 
auroral  oval.  However,  in  assuming  symmetry  between  the 
northern  and  southern  hemispheres,  they  implicitly  solved 
(7)  only  for  equinox  cases.  In  order  to  use  their  formulation 
with  as  little  modification  as  possible,  a  geomagnetic  storm 
near  the  equinox  was  chosen  to  illustrate  our  computations. 

[1 8]  Figure  2  provides  two  examples  of  N-C  model 
solutions  for  potential  distributions  along  the  dusk-dawn 
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Figure  3.  Electric  potentials  from  60°  magnetic  latitude  to 
the  north  pole  with  the  J1  current  shown  as  a  pair  of  gray 
crescents.  The  magnitude  of  J1  is  determined  by  the  Hill- 
Siscoe-Ober  model  and  conditions  in  the  interplanetary 
medium.  In  this  case,  there  is  no  region-2  current  system. 
Contour  values  are  in  kilovolts. 


meridian.  These  plots  were  obtained  using  the  ionospheric 
conductance  as  given  in  (16)  below.  With  region- 1  and 
region-2  current  distributions  described  above,  we  set  J1  = 
1.0  MA  and  J2  =  0.5  MA  (open  circles)  or  0.0  MA  (open 
triangles).  The  plots  in  Figure  2  show  that  with  J2  =  0.5  MA, 
$pc(N-C)  =  107  kV,  and  with  J2  =  0  MA,  $‘PC(N-C)  = 
153  kV.  In  the  present  configuration,  we  are  simply  looking 
at  the  response  of  a  two-dimensional  conductive  ionosphere 
to  freely  chosen  J1  and  J2  FAC.  In  section  5  a  transpolar 
potential  and  a  ring  current  is  added  to  the  model. 


Figure  4.  The  electric  potential  from  the  magnetic  equator 
to  the  north  magnetic  pole  computed  from  the  Nopper- 
Carovillano  model  for  the  condition  described  in  the 
Figure  3.  This  highlights  the  direct  penetration  of  the  polar 
electric  fields  to  equatorial  latitudes. 


[19]  A  typical  ionospheric  N-C  polar  potential  pattern  is 
given  in  Figure  3  together  with  a  pair  of  shaded  regions 
showing  where  the  J1  current  enters  and  leaves  the  iono¬ 
sphere.  The  equipotential  lines  bend  as  they  cross  the  polar 
region  from  the  noon  to  midnight  sectors  due  to  the  gradient 
in  the  conductance.  The  model  has  not  been  adapted  to 
show  the  influence  of  IMF  By  [Weimer,  1995]  which  was 
relatively  small  between  1700  and  2400  UT  on  6  April 
2000. 

[20]  For  weak  J2  the  electric  potentials  extend  equator- 
ward  of  the  auroral  boundary  (\m~63°).  Figure  4  shows 
that  the  ionospheric  potential  driven  by  the  interplanetary 
medium  extends  to  the  magnetic  equator.  Because  the 
conductance  in  the  two  poles  is  assumed  to  be  the  same, 
contributions  from  the  two  hemispheres  are  nearly  the  same 
and  the  N-C  electric  field  at  the  equator  has  only  an 
azimuthal  component.  With  Nopper  and  Carovillano 
[1978],  we  note  that  Wolf  [1910]  demonstrated  the  necessity 
of  an  eastward  electric  field  component  in  the  early  evening 
sector  to  maintain  current  continuity  across  the  conductance 
gradient  at  the  dusk  meridian.  In  Appendix  A,  we  illustrate 
some  ionospheric  features  obtainable  with  N-C. 

5.  Siscoe  Ring  Current  Model 

[21]  Since  one  goal  is  to  determine  the  time-dependent 
penetration  electric  fields  at  lower  latitudes,  we  need  a 
model  for  J2.  For  present  purposes  we  use  that  of  Siscoe 
[1982]  (S-RC),  as  seen  in  Figure  5.  Here  $PC  drives  a 
current  J1  that  splits  into  a  transpolar  component  and  J2' 
that  partially  closes  across  Rs  at  lower  latitudes  and  partially 
across  an  inductor-equivalent  L  in  the  magnetosphere, 
which  represents  the  ring  current.  To  get  a  better  handle 


L 


Figure  5.  Taken  from  Siscoe  [1982],  Rp  is  the  average 
resistance  across  the  polar  cap.  B  and  C  locate  the  upward 
and  downward  legs  of  Jl,  respectively.  J2  enters  the 
ionosphere  at  point  A  and  exits  the  ionosphere  at  D.  The 
ionospheric  resistance  between  A  and  B  and  C  and  D  is  RA. 
Rs  represents  the  effective  resistance  of  the  ionosphere  at 
lower  latitudes  to  closing  current  J2'-J2.  Therefore  in  this 
model  Jl  also  provides  sufficient  current  to  close  at  lower 
latitudes,  across  the  poles,  and  J2  into  the  magnetosphere. 
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on  L,  let  us  visualize  a  fraction  of  $PC  being  applied  across 
closed  field  lines  in  the  magnetotail  corresponding  to  points 
A  and  D  in  Figure  5.  As  is  well-known,  this  results  in  an 
earthward  E  x  B  drift  with  the  ions  gradient-curvature 
drifting  westward.  Hines  [1963]  showed  that  gradient  drift 
causes  perpendicular  heating  of  westward  drifting  ions  at  a 
rate  qVd  •  E,  where  q  is  the  ion  electrical  charge  and  Vd  the 
ion  drift  velocity  due  to  the  magnetic  field  gradient.  The 
ring  current  therefore  has  inductor-like  properties,  storing 
energy  in  the  ion  gyromotion  instead  of  the  magnetic  field. 
The  parameter  L  in  the  S-RC  model  (Figure  5)  reflects  this 
property.  The  Siscoe  [1982]  model  is  intended  to  describe 
the  main  phase  of  a  magnetic  storm. 

[22]  In  Figure  5,  Kirchhoff’s  law  is  satisfied  in  each  of  the 
three  loops,  I,  II,  and  III,  as  now  shown. 

I.  (Jl  -  J2')Rp  =  $pC,  J1  =$PC/Rp  +  J2'  (8) 


II.  (J2'  -  J2)Rs  +  2RaJ2'  =  $pc  (9) 


III.  LdJ2/dt  +  (J2  -  J2')Rs  =  0  (10) 

Solve  (9)  for  J2',  find  J2  -  J2',  and  then  substitute  into  (10). 

LdJ2/dt  +  2RaRsJ2/(Rs  +  2Ra)  =Rs  $/,c/(Rs  +  2RA)  (11) 
Equation  (11)  is  easily  solved. 

J2  =  *«7(2Ra)(1  -  exP( — t/x))  (12) 


t  =  L(Rs  +  2RA)/(2RARs)  (13) 

for  a  step  like  cj>pc  -  J2/  is  found  as 

J2'  =  9kA Rs  +  2Ra)[1  +  Rs/(2RA)(1  -  exp(— t/r))].  (14) 

Therefore  from  (8) 

J1  =  $pc  [1/Rp  +  l/(Rs  +  2Ra)  +  Rs/(2RA(Rs  +  2RA)) 

•  (I  —  exp(— t/r))]  (15) 

[23]  The  first  term  on  the  right-hand  side  of  (15)  denotes 
the  part  of  J 1  crossing  the  polar  cap,  and  the  next  two  terms 
represent  the  part  of  J1  flowing  toward  lower  latitudes 
including  J2  (i.e.,  the  fraction  of  J2’  that  is  diverted  to  the 
magnetosphere).  In  the  limit  of  large  t,  J1  =  $PC/Rp  +  J2. 
The  transpolar  potential  drives  whatever  currents  are  re¬ 
quired  (region- 1)  into  the  dawn  ionosphere  and  out  of  the 
dusk  ionosphere.  Part  of  this  current  directly  crosses  the 
polar  cap  and  the  rest  flows  equatorward,  splitting  into 
region-2  current  and  current  that  flows  through  the  low- 
latitude  and  midlatitude  ionosphere.  If  the  region-2  current 
system  is  fully  developed,  it  tends  to  shield  (cancel)  the 
region- 1  electric  field  and  currents,  especially  at  lower 
latitudes.  This  model  is  consistent  with  satellite  magnetom¬ 
eter  measurements  [Zmuda  and  Armstrong,  1974a,  1974b; 
Iijima  and  Potemra,  1976;  Smiddy  et  ah,  1980]  that  see 
region- 1  and  region-2  currents  as  adjacent  pairs,  a  down- 


Table  1.  Transpolar  Potentials  and  Equatorial  Electric  Fields 
From  the  N-C  Model,  ct  =  0,  P  =  3.8” 


£pp,  S 

Ea.  s 

£ 

E£p,p,  S 

Qpc  b  kV 

Eeql, 

mV/m 

Eeq2, 

mV/m 

4.8 

0 

Jl 

2.01 

=  1  MA 

9.8 

J2  =  0 

101.3 

1.2 

-2.7 

5.6 

3 

2.44 

13.7 

72.9 

1.0 

-1.9 

5.9 

4 

2.51 

15.0 

71.9 

0.8 

-1.6 

7.0 

8 

2.86 

20.0 

50.0 

0.4 

-0.8 

9.1 

16 

3.43 

31.2 

32.0 

0.2 

-0.5 

“With  the  average  polar  cap  Pedersen  conductance  Ep  defined  as  a 
function  of  the  auroral  conductance  SA,  the  transpolar  potential  $PC  and 
the  corresponding  equatorial  electric  fields  Eeql  and  Eeq2  are  found  for 
J1  and  J2.  Eeql  and  Eeq2  are  taken  as  the  best  estimates  of  the  pre¬ 
reversal  maximum  and  minimum,  respectively.  Eeq2  is  found  with  J1  =  0, 
J2  =  1  MA. 


ward  J1  and  an  upward  J2  on  the  dawnside,  and  an  upward 
Jl,  downward  J2  on  the  duskside. 

6.  Combined  Model 

[24]  Using  N-C,  it  is  possible  to  reduce  the  distributed 
polar  current-voltage  relation  to  a  “single-wire”  with  an 
effective  conductance  between  the  maxima  and  minima 
voltages  as  shown  in  Figure  2  and  as  described  above.  This 
allows  us  to  combine  the  N-C  and  H-S-0  models  as  the 
corresponding  voltages  are  applied  across  the  same  “wire.” 
The  N-C  model  is  used  to  determine  E,Ep,  as  stated,  by 
finding  the  ratio  ]\/A>PC(N-C),  usually  with  Jl  =  1  MA  and 
J2  =  0  MA.  The  parameter  6,  partially  depends  on  the 
geometry  of  the  inputted  currents  and  partially  on  the 
ionosphere  presenting  an  impedance  equivalence  of  two 
parallel  resistors.  The  average  Pedersen  conductance  E£ 
can  be  independently  found  by  a  Monte  Carlo  averaging 
of  the  model  conductance  (17)  over  the  polar  cap.  By  using 
this  approach  we  find  £  =  2.0-3.4  (see  Table  1).  Again, 
consistency  requires  that  the  same  value  of  £,e£  be  used  in 
determining  4vc(H-S-0)  from  (2). 

[25]  The  ionospheric  conductance  model  can  be  made 
more  realistic  by  incorporating  its  global  dependence  on 
the  EUV  flux  (via  the  proxy  of  the  F10.7  solar  radio  flux). 
The  global  Pedersen  conductance  Ep(0,  ip)  as  given  by 
Nopper  and  Carovillano  [1978]  is  shown  in  (16). 

Ep(9,  ip)  =  0.3  +  5.0exp  j— [cos-1  (sin(0)cos(ip))]'/1.8O4| 

+  auroral  terms  (16) 

auroral  terms  =  EA  exp(-  (0  -  0a)2/cta);  Ea  =  3  S,  aA  = 
0.57°,  0A  =  20°.  This  expression  can  be  extended  to  arbitrary 
FI 0.7  values  by  using  the  parameters  a  and  |3. 

Ep(9,  ip)  =  Ep^a  +  (3exp|— [cos-1  (sin(9)  cos(ip))]2/1.804|^ 

+  auroral  terms  (17) 

The  average  polar  conductance  Ep  due  to  solar  EUV 
(F10.7)  is  defined  in  (3).  This  implies  that  for  the  6-7  April 
2000  time  period  when  F10.7  =  177,  Ep  =  4.9  S.  We  now 
adjust  a  and  (3,  as  given  in  (17),  to  obtain  an  appropriate 
conductance  value  at  the  subsolar  point.  For  example,  if  a  = 
0  and  P  =  3.8,  one  obtains  a  subsolar  conductance  of  19  S 
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[Senior  and  Blanc,  1984].  The  a  parameter  is  set  to  zero  to 
emphasize  the  local  time  dependence  of  (17).  Consistency 
requires  Ep  =  (Ep(0,  ip)),  where  the  brackets  denote 
averaging  over  the  polar  cap.  We  find  that  (Ep(0,  ip))  = 
4.8  S  in  comparison  with  the  value  of  4.9  S  obtained  from 
(3).  It  is  concluded  that  with  these  values  for  a  and  0,  (17)  is 
a  reasonable  representation  of  the  global  ionospheric 
conductance  on  6-7  April  2000.  Having  defined  the  global 
ionospheric  conductance  we  can  directly  relate  the  iono¬ 
spheric  electric  fields  at  any  location  in  the  ionosphere  to 
Esw.  In  Appendix  A  we  show  some  of  the  high-  and  low- 
latitude  features  of  N-C. 

7.  Application  of  the  Combined  Model 

[26]  In  the  present  approach  there  are  three  models  that 
form  a  basic  solar  wind-magnetosphere-ionosphere  model. 
According  to  the  Hill-Siscoe  (H-S)  model,  the  solar  wind  at 
LI  detennines  the  transpolar  potential  <bPC  over  the  polar 
caps,  consistent  with  the  polar  conductance  predicted  by  the 
ionospheric  model  (N-C).  The  results  are  shown  to  be  in 
good  agreement  with  both  the  empirical  Weimer  [2005] 
transpolar  potentials  and  DMSP  measurements  for  the  6-7 
April  2000  magnetic  storm.  The  region- 1  and  region-2  (J2) 
currents  are  estimated  from  the  Siscoe  (S-RC)  ring-current 
circuit  model  which  is  driven  by  §PC.  The  resistor  values  in 
S-RC  are  determined  by  N-C  which,  when  combined  with 
the  global  potential  solution,  make  it  possible  to  estimate 
the  time  profile  of  the  equatorial  penetration  electric  field 
during  the  storm’s  main  phase.  Applying  this  approach  to 
the  6-7  April  2000  storm,  we  find  that  shielding  occurs 
within  1  hour  of  onset. 

[27]  The  J2  rise  time  is  estimated  using  the  Dessler- 
Parker-Scopke  (DPS)  [Dessler  and  Parker,  1959;  Scopke, 
1966]  relation.  In  this  manner,  expressions  for  J1  and  J2  are 
found  that,  when  combined  with  the  equatorial  electric  fields 
estimated  from  N-C,  give  a  time-dependent  expression  for 
shielding.  It  is  to  be  noted  that  the  obtained  shielding  times 
tend  to  be  much  shorter  than  the  J2  rise  times. 

7.1.  High  Latitudes 

[28]  A  key  feature  therefore  is  to  accurately  define  the 
polar  ionospheric  conductance  during  the  6-7  April  2000 
magnetic  storm.  We  do  this  by  considering  five  values  of 
auroral  conductance,  EA  =  0,  3,  4,  8,  and  16  S,  in 
determining  fEp.  The  five  corresponding  values  of  6,Ep, 
when  put  into  (2),  give  five  curves  for  TVfjH-S-O)  (see 
Figure  6c)  throughout  the  6-7  April  2000  magnetic  storm. 
These  curves  are  then  compared  with  the  experimental 
$/>c(DMSP)  values  for  the  transpolar  potential  as  shown 
in  Figure  6B.  We  then  choose  the  value  of  EA  for  which 
<3?/jc(H-S-0)  and  $/>c(DMSP)  are  the  best  match.  In  this 
case  EA  =  3  S.  This  value  of  EA  is  also  consistent  with  the 
observed  electron  precipitation  energy  using  the  formula  of 
Robinson  et  al.  [1987].  The  “standard”  value  for  (.E^ 
during  the  6-7  April  2000  storm  is  therefore  13.7  S  as 
given  in  Table  1. 

[29]  By  using  N-C  to  define  the  conductance  (;Ep  for  the 
H-S-0  model  (2),  we  find  good  agreement  with  the  Weimer 
[2005]  empirical  potential  model  for  this  magnetic  storm,  as 
shown  in  Figure  6a.  Note  that  the  Weimer  model  is  based  on 
DE-2  satellite  data  taken  between  August  1981  and  March 


1983.  During  this  period  the  value  of  F10.7  was  178  ±  33 
compared  with  177  for  the  present  magnetic  storm  which 
indicates  the  comparison  is  being  made  under  similar 
conditions  of  global  EUV  conductance.  In  the  Weimer 
[2005]  model  a  parameterized  exponential  function  is  used 
to  model  saturation  effects  in  contrast  to  the  H-S-0  expres¬ 
sion  given  in  (2).  Note  that  the  agreement  between  the 
Weimer  and  H-S-0  models  is  dependent  on  £E£,  as  calcu¬ 
lated  by  N-C.  Therefore  Figure  6a  represents  an  overall 
consistency  between  these  three  models. 

[30]  Experimental  values  of  the  transpolar  potentials 
shown  in  Figure  6b  were  derived  from  data  taken  on  board 
DMSP  satellite  FI 3.  During  the  6-7  April  2000  magnetic 
storm  satellite  F 1 3  approached  the  center  of  the  N-C  potential 
cells  shown  in  Figure  3,  particularly  during  northern  passes. 
Under  these  conditions  good  agreement  was  found  between 
the  measured  potentials  and  the  values  predicted  by  N-C. 
Therefore  <F/jc(DMSP)  Was  estimated  from  the  potential 
minima  and  maxima  [Rich  and  Hairston,  1994]  and  plotted 
as  open  squares  in  Figure  6b.  Note  that  the  DMSP  results 
generally  agree  with  those  of  Weimer  in  Figure  6a. 

[31]  Figure  6c  shows  five  time  histories  of  4Vc(H-S-0)  as 
determined  by  the  ACE  satellite  solar  wind  measurements. 
The  five  curves  represent  the  five  values  of  auroral  con¬ 
ductance  Ea  as  given  in  Table  1  with  the  solar  wind 
parameters  shifted  to  later  times  consistent  with  the  solar 
wind  velocity.  The  five  curves  in  Figure  6c  represent  the 
effect  of  enhanced  auroral  conductance  on  $PC(H-S-0).  In 
other  words,  auroral  precipitation  suppresses  $pC(H-S-0). 
Roughly,  there  is  about  a  16  kV  decrease  in  <FpC(H-S-0)  for 
each  1  S  increase  in  Ea  between  0  and  3  S. 

7.2.  Equatorial  Regions 

[32]  Now  let  us  derive  the  equatorial  electric  field  Eeq 
using  the  N-C  model.  The  N-C  model  is  treated  as  a 
numerical  Green’s  ftinction.  That  is,  for  a  fixed  geometry 
and  global  conductance  the  potential  and  electric  fields  at 
specific  locations  are  determined  by  independently  input¬ 
ting  1  M A  values  for  J1  and  J2.  The  corresponding  potential 
and  electric  field  values  at  these  locations  for  arbitrary  J1 
and  J2  may  then  be  found  by  scaling  and  superposition.  We 
denote  Eeql  for  the  prereversal  enhancement  electric  field 
from  J1  and  Eeq2  as  the  corresponding  minimum  (negative) 
electric  field  from  J2. 

[33]  Table  1  lists  the  results  for  <S>pc(N-C)  and  the 
equatorial  electric  fields  Eeql  and  Eeq2  using  the  five 
values  of  the  auroral  conductance  EA,  0,3,  4,  8,  and  16  S 
[Gonzales  et  al.,  1983],  For  each  case  a  Jl  =  1  MA  current 
is  put  into  the  N-C  model  as  being  evenly  distributed  into 
and  out  of  two  120°  crescents  at  \m  =  12°  centered  at  the 
dusk  and  dawn  meridians.  Consistent  with  the  standard 
convention  J1  enters  the  ionosphere  on  the  momingside 
and  exits  on  the  duskside.  The  crescents  have  a  width  of  3° 
so  that  a  1  MA  current  corresponds  to  a  parallel  current 
density  of  0.72  qA/m2.  The  region  2  current  J2  is  placed  at 
66°,  just  equatorward  of  Jl,  also  with  a  width  of  3°  and 
extending  120°  in  longitude.  J2,  of  course,  has  the  opposite 
sense  of  entering  and  exiting  the  ionosphere  as  does  Jl.  In 
the  present  instance  J2  =  1  MA  corresponds  to  a  parallel 
current  density  of  0.56  pA/m2. 

[34]  It  should  be  pointed  out  that  although  for  modeling 
purposes  we  keep  the  latitudes  of  J 1  and  J2  fixed  throughout 
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Figure  6.  (a)  Comparison  of  the  transpolar  potential  $PC  from  the  Weimer  [2005]  model  for  the 
6-7  April  2000  magnetic  storm  with  that  from  the  H-S-0  model  for  zero  and  three  Siemens  auroral 
conductance.  The  inserted  plot  ranges  from  1600  UT  on  6  April  to  0100  UT  on  7  April  (in  one  hour 
intervals)  in  order  to  highlight  peak  details.  The  y-axis  ranges  from  150  to  250  kV  in  25  kV  steps,  (b)  The 
EA  =  3  S  H-S-0  curve  is  compared  with  corrected  DMSP  measurements  of  the  transpolar  potential 
during  the  same  storm  using  the  N-C  model.  We  use  only  northern  passes  of  satellite  FI 3  as  these  go 
more  directly  through  the  potential  cell  centers,  (c)  The  H-S-0  transpolar  potential  as  calculated  at 
various  levels  of  auroral  conductance  as  shown.  Note  the  sensitivity  of  the  results  at  lower  values  of  SA. 
On  the  basis  of  the  electron  precipitation  energy,  EA  =  3  S  is  considered  to  be  the  correct  theoretical 
curve. 


the  stonn,  it  is  recognized  that  during  the  course  of  the 
storm  the  auroral  zone  moved  to  lower  latitudes.  The  auroral 
oval  boundaries  expand  equatorward  after  significant  net 
open  flux  has  been  added  to  the  polar  cap  [SAcoe  and 
Huang,  1985]  and  the  plasma  sheet  convects  close  to  Earth 
[Croaker  and  Siscoe,  1983],  This  expansion  occurred  with  a 
lag  of  about  3  hours  between  the  beginning  of  the  storm  and 
the  time  of  significant  equatorward  movement  [Huang  and 
Burke,  2004]. 

[35]  From  Table  1  and  Figure  6c  it  is  also  seen  that  a 
higher  auroral  conductance  depresses  §>PC  $PC(N-C)  at 
constant  Jl.  Put  another  way,  the  polar  ionosphere  with  a 
higher  conductance  demands  a  larger  current  in  order  to 
sustain  the  same  value  of  <&PC.  This  causes  saturation  of  Tpc 
to  occur  at  lower  values  as  can  be  seen  from  (2). 

[36]  We  now  turn  to  the  penetration  electric  field  as 
calculated  by  N-C.  The  conductance  gradients  arising  from 
(17)  cause  Jl  to  create  a  “prereversal”  enhancement  in  Eeq 
on  the  eveningside,  although  significantly  displaced  from  its 


proper  location  at  1 800  MLT  [Fejer,  1997].  It  is  found  that 
this  displacement  is  due  to  the  approximate  nature  of  (17), 
which  does  not  contain  the  small-scale  features  of  the 
terminator  conductance.  The  equatorial  electric  fields  are 
essentially  fringe  electric  fields  and  therefore  are  sensitive 
to  source  location.  They  are  also  difficult  to  calculate 
accurately  because  of  strong  conductance  gradients  near 
the  equator  from  dip  angle  effects  and  density  gradients. 
However,  one  modification  including  density  gradient 
effects  near  the  tenninator  was  used  with  the  results  shown 
in  Figure  7  for  Eeql.  There  is  a  similar  curve  for  Eeq2, 
but  with  the  opposite  sign  and  larger  amplitude.  For 
modeling  purposes  we  use  the  extremum  (prereversal) 
values  of  Eeql  and  Eeq2  as  found  from  N-C  with  Jl 
and  J2  separately  set  equal  to  1  MA,  as  shown  in  Table  1 . 
One  sees  from  Table  1  that  Eeq2  ~  -2  Eeql,  which 
reflects  sensitivity  to  the  location  of  the  source  currents, 
since  J2  is  closer  to  the  equator.  It  is  well  known  that 
during  a  magnetic  storm  both  J 1  and  J2  move  closer  to  the 
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Figure  7.  The  equatorial  electric  field  Eeql  as  a  function 
of  MLT,  J1  =  1  MA,  Xm  =  72°,  and  J2  =  0  MA.  This 
calculation  included  an  additional  density  gradient  at  the 
tenninator  which  placed  the  pre-reversal  enhancement  near 
1 800  MLT.  The  calculation  was  performed  using  a  banded- 
matrix  method  as  described  by  Press  et  al.  [1992],  A  similar 
curve  is  found  for  Eeq2  (J1  =  0  MA,  J2  =  1  MA,  Xm  =  66°), 
but  of  the  opposite  sign  and  with  a  larger  amplitude.  Note 
that  the  integral  of  Eeq  around  the  equator  is  zero,  since  curl 
E  =  0. 


equator.  However,  it  is  found  that,  at  present,  because  the 
net  Eeq  is  a  small  difference  of  two  relatively  large 
numbers  the  error  in  the  penetration  electric  field  Eeq 
becomes  too  large  if  J 1  and  J2  move  substantially  equator- 
ward.  This  feature  will  be  addressed  in  future  work. 

[37]  A  decrease  in  magnitude  of  both  Eeql  and  Eeq2 
occurs  as  the  auroral  conductance  £A  increases  (Table  1). 
On  the  other  hand,  both  Eeq  1  and  Eq2  are  linearly  related  to 
J1  and  J2,  respectively,  as  expected  from  Poisson’s  equa¬ 
tion.  Therefore  one  can  apply  superposition  in  finding  the 
net  Eeq  for  fixed  EA.  Note  also,  however,  that  the  equatorial 
ionosphere  is  three-dimensional,  not  two-dimensional  as 
treated  here.  This  feature  may  make  a  significant  change 
in  the  present  values  of  Eeq. 

[38]  By  using  the  results  from  Table  1  and  superposition 
we  have  the  following  relation  for  each  specified  value  of 
Ea. 

Eeq(mV/m)  =  Eeql  Jl(MA)  +  Eeq2  J2(MA)  (18) 

Equation  (18)  gives  a  simple  representation  of  the 
“prereversal”  electric  field  as  a  function  J1  and  J2. 
However,  as  shown  above,  S-RC  indicates  that  J1  and  J2 
are  coupled,  so  we  have  the  more  complex  expression  (19) 
for  the  penetration  electric  field  Eeq(t)  from  inputting  (12) 
and  (15)  into  (18). 

Eeq(t)  =  $/>c(l/Rp  +  l/(Rs  +  2Ra)  +  Rs/(2RA(Rs  +  2RA)) 

•  (1  -  exp(— t/T)))Eeql  +  $>pC/(2RA) 

•  (1  -  exp(-t/i-))Eeq2  (19) 


7.3.  Using  N-C  to  Define  the  Resistors  in  S-RC 

[39]  Owing  to  the  tensor  nature  of  the  ionospheric  conduc¬ 
tance  (6)  finding  values  for  the  resistors  shown  in  Figure  5 
is  not  obvious.  However,  since  only  the  diagonal  compo¬ 
nents  of  the  conductance  tensor  contribute  to  energy  dissi¬ 
pation  (J-E  ^  0),  the  reciprocals  of  (See)  and  (Sw)  can  be 
used  to  determine  Rp,  RA,  and  Rs,  where  the  bracket 
denotes  an  average  value.  From  Figure  5  it  is  immediately 
seen  that  an  applied  voltage  across  BC,  which  represents  the 
polar  cap,  is  equivalent  to  placing  a  voltage  across  two 
resistors  in  parallel,  with  L  removed  from  the  circuit.  For 
two  resistors  in  parallel  the  total  conductance  is  the  sum  of 
the  two  individual  conductances.  This  is  partially  why  the  £, 
parameter,  in  the  present  treatment,  varies  between  2  and 
3.4  in  Table  1. 

[40]  The  centers  of  the  polar  potential  cells  coincide  with 
the  location  of  the  region-1  (Jl)  currents.  Therefore  points  B 
and  C  in  S-RC  (Figure  5),  in  the  present  treatment,  corre¬ 
spond  to  72°  magnetic  latitude  in  the  N-C  model.  In  a 
similar  manner,  points  A  and  D  (J2)  in  Figure  5  correspond 
to  66°  in  magnetic  latitude.  The  cross-cap  Pedersen  current 
flows  mainly  in  the  0  direction  so  that,  using  (17)  and  (6)  for 
magnetic  latitudes  between  72°  and  90°,  Rp  =  (See)-1  = 
0.174  ft.  Throughout,  SA  is  set  equal  to  3  S. 

[41]  We  now  want  to  estimate  2RA  in  Figure  5.  The 
average  conductance  (See)  between  66°  and  90°  in 
magnetic  latitude  (points  A  and  D)  is  found  such  that 
Rp  +  2Ra  =  (See)-1  =  0.210  ft.  A  value  of  0.174  ft  has 
just  been  found  for  Rp  so  that  2RA  =  0.036  ft. 

[42]  To  determine  Rs,  we  consider  a  circuit-equivalent 
equatorward  of  66°.  The  low-latitude  currents  close  by 
flowing  equatorward  from  66°  in  latitude  on  the  dawnside 
and  then  azimuthally  in  two  paths,  pass  noon  and  midnight, 
where  they  flow  poleward  to  close  on  the  duskside.  The 
S-RC  resistor  Rs  therefore  is  composed  of  two  parallel 
azimuthal  resistors,  one  for  the  dayside  and  one  for  the 
nightside,  and  two  meridian  resistors  representing  dusk  and 
dawn.  We  find  from  the  N-C  model  (Se0)  (dawn)  =  30.65  S, 
(See)  (dusk)  =  30.22  S,  (Sw)  (day)  =  16.56  S,  (Sw) 
(night)  =  2.75  S.  As  expected,  most  of  the  current  closes 
on  the  dayside.  The  value  for  the  resistor  Rs  is  now  given 
in  (20). 

Rs=  1/30.65  +  1/30.22  +  1/(16.56  +  2.75)  =  0.1 17  Q 

(20) 

In  summary,  the  following  values  for  the  three  resistors  are 
Rp  =  0.174  ft,  2Ra  =  0.036  ft,  and  Rs  =  0.117  ft. 

7.4.  Using  DPS  to  Estimate  J2  Rise  Time 

[43]  The  rise  time  for  J2  may  be  estimated  from  the 
Dessler-Parker-Scopke  (DPS)  [Dessler  and  Parker,  1959; 
Scopke,  1966]  relation  that  relates  Dst  with  the  total  ring 
current  energy.  Including  Earth  conductivity  effects,  a  1 00  nT 
decrease  in  Dst  corresponds  to  2.8  x  101 5  J  of  ring  current 
energy  [Kivelson  and  Russell,  1997].  In  tenns  of  the  present 
model  this  relation  can  be  recast  as 

'/2L  J22  =2.8  x  10i3|Dst|.  (21) 
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Figure  8.  Model  time  history  of  J1  and  J2  currents  during  the  6-7  April,  2000  magnetic  storm  using 
(15)  and  (12),  respectively,  based  on  the  model  described  by  Siscoe  [1982]  and  as  shown  in  Figure  5 
with  resistor  values  Rp  =  0.17  Q,  RA  =  0.036  0,  and  Rs  =  0.11  O.  These  values  were  found  using  N-C. 
We  also  used  TVcfH-S-O)  (2)  with  SA  =  3  S,  as  shown  in  Figure  6b,  in  obtaining  J1  and  J2  from  (15) 
and  (12). 


The  J2  rise  time  t  is  substituted  for  L  in  (21)  using  (13). 

t  =  1.12  x  1 014 |Dst| (Rs  +  2Ra)Ra/(Rs$£c)  (22) 

Obviously,  (22)  is  inconsistent  with  (12),  as  the  latter 
assumes  t  is  constant.  Flowever,  we  will  use  (22)  in  an 
empirical  manner  to  estimate  t  in  order  to  determine  J 1  (t) 
(15)  and  J2(t)  (12).  In  (22)  we  assumed  that  J2  has  reached 
its  maximum  value  of  <$>pC/(2RA)  =  5.5  MA  by  the  time  Dst 
reached  its  minimum  value  (—290  nT).  The  argument  being 
made  that  under  these  conditions  (22)  is  a  reasonable 
approximation  to  the  true  rise  time.  The  following  value  is 
found  for  t. 

t  =  5.5h  L  =  545H  (23) 

The  J2  rise  time  t  may  be  longer  than  the  shielding  time,  as 
the  shielding  time  also  depends  on  $/>c/(2RA),  as  will  be 
shown  below.  Therefore  t  is  somewhat  longer  than  the 
shielding  times  found  by  0.5  hours  [Senior  and  Blanc, 
1984],  ~1  hour  [Fejer,  1997],  or  3-300  min  [Richmond  et 
al.,  2003].  The  total  ring  current  energy  at  Dst  minimum  is 
6.9  x  1 015  J  by  (21)  which  is  within  the  usual  range  of 
1 0 1 5  —  1 0 1 6  J  [Nagatsuma,  2002],  Circuit  analysis  regarding 
magnetospheric-ionospheric  coupling  has  also  been  carried 
out  by  Bostrom  [1964],  Weimer  [1994],  and  Block  et  al. 
[1998], 

7.5.  Shielding 

[44]  In  Figure  8,  J1  (15)  and  J2  (12)  are  plotted  using  the 
resistor  values  found  above  and  the  H-S-0  transpolar  poten¬ 
tial  values  from  Figure  6b.  Both  J1  and  J2  in  Figure  8 


reach  large  values  later  into  the  storm  as  observed  by 
Huang  and  Burke  [2004].  Now  combining  Jl(t)  and  J2(t) 
with  the  values  Eeql  =  1.0  mV/m  and  Eeq2  =  -2.0  mV/m 
from  Table  1,  the  time-dependent  penetration  electric  field 
Eeq(t)  (19),  as  shown  in  Figure  9,  is  found.  The  value  Eeq  = 
1  mV/m  corresponds  to  an  upward  drift  velocity  at  Jica- 
marca  of  40  m/s  [Fejer,  1997],  The  horizontal  line  at  Eeq  = 
1.25  mV/m  (Figure  9)  then  corresponds  to  an  upward 
velocity  of  ~50  m/s,  a  value  found  by  Fejer  [Fejer  et  al., 
1999]  and  Whalen  [Whalen,  2001]  as  the  threshold  for  the 
generation  of  strong  spread  F  and  the  onset  of  bubbles.  Note 
from  Figure  8  that  at  the  beginning  of  the  storm  (1643  UT) 
J1  increases  much  faster  than  J2,  resulting  in  an  initial 
positive  spike  in  Eeq(t)  (Figure  9).  Five  bubbles  were 
observed  on  the  DMSP  satellites  between  1818  and 
2037  UT  on  06  April  2000.  These  are  indicated  by  symbols 
at  the  bottom  of  Figure  9  and  are  explicitly  displayed  in 
Figure  10.  Looking  at  Figure  9,  it  is  seen  that  some  of  these 
earlier  bubbles  could  be  related  to  the  initial  ramping  of  Eeq. 
Later  into  the  storm  other  causes,  such  the  neutral  wind 
dynamo  are  probably  responsible  for  bubbles.  The  timing  of 
Eeq  dropping  below  the  Fejer-Whalen  threshold  is  very 
dependent  on  the  equivalent  resistor  values  found  for  S- 
RC.  For  example,  if  we  increase  RA  by  a  factor  of  two,  then 
Eeq  drops  below  the  Fejer-Whalen  threshold  2.4  hours  later, 
as  seen  in  Figure  9.  The  increase  in  shielding  time  is  due  to 
the  effect  of  RA  on  J2  as  is  easily  determined  from  (12) 
and  (22).  Later  in  the  storm,  magnetometers  on  board  the 
DMSP  satellites  observe  Region  2  precipitation  currents  at 
low  magnetic  latitudes  ~50°.  These  currents  strongly 
imply  shielding  consistent  with  Figures  8  and  9.  Therefore 
the  later  bubbles  are  not  produced  by  penetration  electric 
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Figure  9.  Time-history  on  6  April  2000  of  Eeq  (~1800  MLT)  based  on  the  values  of  Jl,  J2  as  given  in 
Figure  8  and  the  values  of  Eeql  and  Eeq2  as  given  in  Table  1  of  text  for  SA  =  3  S.  The  solid  line 
represents  our  best  estimate.  An  increase  in  RA  suppresses  shielding,  as  does  a  decrease  in  Rs,  the  low- 
latitude  resistance.  All  resistance  values  are  in  Ohms.  In  terms  of  Figure  5  an  increase  in  RA  lengthens  the 
J2  rise-time  as  well  depresses  the  J2  amplitude  (12).  A  decrease  in  Rs  diverts  current  through  the  low- 
latitudes  that  otherwise  would  have  contributed  to  J2.  The  curves  are  continued  beyond  24:00  to  show  the 
relation  with  the  two  later  bubbles.  The  solid  line  represents  the  Fejer- Whalen  [Fejer  et  al,  1999; 
Whalen ,  2001]  threshold  for  the  onset  of  strong  spread  F  (i.e.,  bubbles). 


fields.  Two  bubbles  were  observed  late  in  the  storm  on  7 
April  at  0058  UT  and  at  0139  UT.  (Given  the  closeness  in 
geographic  longitude  of  the  two  bubbles,  they  are  probably 
the  same  bubble  measured  by  two  different  satellites.) 
Mechanisms,  such  as  the  ionospheric  winds  discussed  by 
Richmond  et  al.  [2003],  may  play  a  key  role  here  [see 
Richmond  et  al.,  2003,  Figure  4].  It  should  be  emphasized 
that  DMSP  spacecraft  do  not  detect  every  bubble  due  to 
their  altitude  (835  km)  and  orbital  period  (101  min). 

8.  Summary 

[45]  Three  models  have  been  combined  to  form  a  solar 
wind-magnetosphere-ionosphere  model.  Our  approach  has 
been  to  use  the  strengths  of  one  model  to  offset  deficiencies 
in  the  other  two.  For  example,  the  ionospheric  model  N-C 
[Nopper,  1978;  Nopper  and  Carovillano,  1978,  1979], 
given  Jl  and  J2,  replicates  the  expected  ionospheric  features 
quite  well.  The  Jl  current  is  driven  by  the  solar  wind  as 
described  by  the  Hill-Siscoe  model  [ Hill  et  al,  1976;  Siscoe 
et  al.,  2002].  The  J2  current,  according  to  the  S-RC  [Siscoe, 
1982]  model,  is  driven  by  Jl  so  that  the  region- 1  and 
region-2  currents  are  coupled  [Siscoe,  1982].  Knowledge 
of  J2  is  necessary  in  order  to  estimate  the  penetration 
electric  field.  From  another  point  of  view,  Hill-Siscoe 
requires  knowledge  of  the  polar  cap  Pedersen  conductance 
to  determine  the  transpolar  potential.  This  is  supplied  by 
N-C.  In  addition,  S-RC  [Siscoe,  1982]  uses  Hill-Siscoe  to 
define  a  transpolar  potential  that  drives  a  simple  LR  circuit, 
the  resistive  element  values  are  defined  from  N-C.  The 


result  is  a  unified  model  that  gives  the  time-dependent 
behavior  of  the  penetration  electric  field  based  on  the  solar 
wind  at  LI. 

[46]  In  more  detail,  we  have  explored  the  applicability  of 
Poisson’s  equation  to  a  two-dimensional  conducting  iono¬ 
sphere  including  dip-angle  effects.  Following  Nopper  and 
Carovillano  [1978],  it  was  found  that  such  a  simple  treat¬ 
ment  yields  two-cell  polar  cap  potentials,  electrojets  (see 
Appendix  A),  as  well  as  giving  estimates  for  the  equatorial 
electric  fields  in  terms  of  region- 1  and  -2  currents.  The  latter 
was  simplified  by  using  the  concept  of  Green’s  functions. 
The  N-C  model  provides  a  systematic  procedure  for  finding 
the  conductance  appropriate  for  the  H-S-0  model  consistent 
with  conductance  values  at  lower  latitudes.  The  resulting 
values  for  the  H-S-0  transpolar  potential  were  found  to  be 
in  very  good  agreement  with  those  of  the  Weimer  [2005] 
empirical  model  (Figure  6a)  and  DMSP  (Figure  6b).  The 
Siscoe  [1982]  ring  current  model  was  used  to  attach  a 
simple  magnetospheric  circuit  to  the  N-C  model.  Conduc¬ 
tance  properties  of  the  N-C  model  were  used  to  define 
resistor  values  in  the  Siscoe  [1982]  model.  In  this  manner, 
one  can  estimate  ring  current  (J2)  response  to  solar  wind 
parameter  changes.  An  equivalent  inductance  L  for  the  ring 
current  was  justified  based  on  the  work  of  Hines  [1963], 
who  showed  that  electrical  energy  was  directly  stored  in  ion 
gyromotion.  However,  while  the  Dessler-Parker-Scopke 
relation  allows  reasonable  J2  rise  times  to  be  estimated, 
there  are  serious  consistency  problems  regarding  L  that 
need  to  be  resolved.  Finally,  in  Figure  9  we  estimated  the 
penetration  electric  field  as  a  function  of  time.  The  solid  line 
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Figure  10.  Observations  of  total  ion  density  at  845  km  by  DMSP  spacecraft  on  6-7  April  2000  as  a 
function  of  the  corrected  geomagnetic  latitude  at  the  subsatellite  point.  The  UT  and  MLT  times  in  the 
upper  right  of  each  frame  and  the  longitude  in  the  lower,  center  are  values  when  the  spacecraft  crossed  the 
magnetic  equator.  All  three  spacecraft  travel  south  to  north  in  the  evening  sector.  The  abscissa  denotes 
magnetic  latitude  with  zero  at  the  center  of  each  plot  and  70°S  and  70°N  denoting  the  limits. 


represents  the  case  for  the  resistor  values  found  from  the 
N-C  model  and  indicates  that  Eeq(t)  stayed  above  the 
Fejer-Whalen  threshold  for  about  an  hour  after  onset.  If, 
however,  one  doubles  the  value  of  RA  in  the  Siscoe  [1982] 
model  from  0.018  S  to  0.036  S,  then  Eeq(t)  is  above  the 
Fejer-Whalen  threshold  for  3.5  hours  as  also  seen  in 
Figure  9.  Therefore  even  in  this  simple  treatment  penetra¬ 
tion  electric  fields  are  dependent  on  nonlocal  effects,  as 
Ra  is  one  measure  of  resistance  at  high  latitudes. 

9.  Discussion 

[47]  Our  coupling  of  three  models,  N-C,  H-S-O,  and 
S-RC  helps  to  explain  the  limited  duration  of  electric-field 
penetration  in  geophysical  terms.  The  sudden  increases  in 
the  pressure  on  the  Earth’s  magnetosphere  and  the  inter¬ 
planetary  electric  field  almost  immediately  drove  a  J1 
current,  as  implied  by  the  H-S-0  model.  The  Jl-flow  occurs 
promptly  because  its  source  is  tied  directly  to  the  magneto¬ 
sheath  generator  via  newly  opened  magnetic  field  lines.  The 
agreement  between  H-S-0  and  Wiemer  [2005],  as  seen  in 
Figure  6a,  using  N-C  under  the  present  assumptions, 


indicates  that  the  Hill-Siscoe  model  accurately  reflects 
the  dependence  of  the  transpolar  potential  on  the  solar 
wind.  Companion  J2  currents  cannot  respond  immediately 
because  their  source  is  distributed  in  the  plasma  sheet  and 
inner  magnetosphere.  The  pressure  gradients  needed  to 
drive  these  currents  [Vasyliunas,  1970]  typically  take 
several  hours  to  grow  [ Garner  et  al.,  2004],  analogous 
to  the  response  of  an  inductor.  At  storm  onset  the  plasma  that 
provides  the  high-pressure  source  of  J2  current  is  spread 
throughout  the  geomagnetic  tail  and  requires  ~3  hours  to 
concentrate  in  the  inner  magnetosphere,  a  timescale  that  is 
consistent  with  the  times  found  above. 

[48]  Bubbles  can  form  in  the  postsunset  equatorial  iono¬ 
sphere  due  to  the  E  x  B  force  pushing  plasma  upwards. 
Whalen  [2001]  showed  that  bubbles  generally  form  when 
the  upward  ion  drift  is  50  m/s  or  greater.  If  the  normal 
diurnal  cycle  of  equatorial  plasma  drifts  [e.g.,  Fejer  et  al., 
1999]  is  augmented  by  electric-field  penetration  from  high 
latitudes  in  sufficient  strength  to  cross  the  50  m/s  (1 .25  mV/m) 
Fejer-Whalen  threshold,  then  EPBs  and  accompanying 
radio-signal  scintillations  are  likely  to  follow  as  shown  in 
Figure  9.  Note,  however,  that  not  all  bubbles  initiated  at 
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lower  altitudes  will  be  observed  on  DMSP  satellites  and  that 
other  causes  such  as  global  winds  driven  by  solar  heating  or 
disturbance  winds  driven  by  high-latitude  heating  and  ion- 
drag  acceleration,  as  mentioned  by  Richmond  et  al.  [2003], 
may  also  cause  bubbles.  A  zonally  propagating  gravity 
wave  can  also  initiate  the  Rayleigh-Taylor  instability  in 
the  bottomside  F  region  [Huang  and  Kelley,  1996a,  1996b], 

[49]  The  6-7  April  2000  storm  started  at  1643  UT 
and  DMSP  spacecraft  observed  EPBs  from  ~  1 8 1 7  UT 
to  ~2040  UT.  Looking  at  Figure  9,  one  sees  that  Eeq 
first  crossed  the  Fejer- Whalen  threshold  at  1655,  less  than 
2  hours  before  the  first  bubble  was  observed  at  1818.  The 
second  set  of  bubbles  were  initiated  at  0058  on  7  April, 
more  than  8  hours  after  Eeq  exceeded  the  Fejer-Whalen 
threshold,  indicating  they  were  probably  caused  by  other 
mechanisms. 

[50]  In  addition  to  the  increasing  strength  of  J2  currents, 
both  J1  and  J2  moved  to  lower  latitudes  between  1900  and 
2200  UT.  More  details  of  this  equatorial  movement  for  the 
6  April  2000  stonn  may  be  seen  in  Figure  4  of  Huang  and 
Burke  [2004],  Early  in  the  storm  (~1707  UT),  DMSP 
observations  near  the  dawn  tenninator  place  J1  at  ~73°. 
As  previously  mentioned,  the  increase  in  both  model  Eeql 
and  Eeq2  as  J 1  and  J2  move  equatorward  compromises  the 
accuracy  of  Eeq,  which  is  dependent  on  the  difference  of 
the  two.  Also,  we  are  presently  improving  the  conductivity 
gradient  near  the  tenninator  to  give  us  more  realistic 
estimates  of  Eeq.  Both  the  location  of  J1  and  the  strength 
of  J2  need  to  be  modeled  as  accurately  as  possible  in 
estimating  the  point(s)  in  a  storm  when  the  upward  plasma 
drifts  reach  the  threshold  for  generating  EPBs. 

[51]  We  have  assumed  that  the  polar  conductance  remains 
constant  throughout  the  stonn  which  allowed  us  to  use  a 
linear  treatment.  However,  this  approximation  may  not  hold 
in  all  cases,  particularly  when  the  electron  precipitation  is 
very  energetic.  The  solution  for  the  general  problem 
involves  an  implicit  relation  for  J1  as  shown  in  (24). 


= _ 3  OP]/,2  +  57.6E.wPff _ 

0.0187£,Ep(Jl)P^6  +  0.036ESW^(J1)  +  P^2 
=  4>pc(H  -  S  -  O)  (24) 

The  dependence  on  J1  enters  through  its  effect  on  auroral 
conductance,  and  therefore  ££p(Jl)  may  be  numerically 
found  by  using  N-C  to  determine  <&PC(N-C)  for  a  range 
ofJl. 

10.  Conclusions 

[52]  1.  There  is  agreement  between  H-S-0  and  the 
Weimer  model  as  seen  in  Figure  6a,  since  F10.7  ~  177 
for  both  the  present  magnetic  stonn  and  during  the  time 
interval  on  which  the  Weimer  model  [Weimer,  2005]  is 
based.  The  similar  values  of  FI 0.7  imply  the  same  EUV 
polar  conductance  was  present  in  both  cases.  The  agreement 
also  indicates  that  the  N-C  model  gives  reasonable  conduc¬ 
tance  values  as  input  to  the  H-S-0  model. 

[53]  2.  The  N-C  and  H-S-0  models  can  be  made  consis¬ 
tent  with  one  another  by  using  N-C  to  define  the  equivalent 


ELECTRIC  FIELD  J1  =  1  MA  J2  =  0  MA  SA=  3  S 


Figure  Al.  This  figure  shows  the  electric  field  at  high 
latitudes  for  the  ionospheric  conductance  used  in  the  text 
(17).  Note  the  positive  charge  on  the  dawnside  and  the 
negative  charge  on  the  duskside. 

“single-wire”  conductance  between  the  potential  extremes 
shown  in  Figure  2. 

[54]  3.  The  N-C  model  may  be  treated  as  a  numeric 
Green’s  function  provided  that  the  ionospheric  conductance 
and  the  J 1 ,  J2  geometries  remain  unchanged.  By  running  the 
two  cases  J1  =  1  MA,  J2  =  0  and  J1  =  0,  J2  =  1  MA  the 
potential  and  electric  field  at  any  point  for  arbitrary  values 
of  J 1 ,  J2  may  be  obtained.  The  N-C  model  can  also  be  made 
more  realistic  by  adding  a  conductance  gradient  at  the 
terminator  based  on  density  data.  This  addition  has  the 
effect  of  shifting  the  prereversal  peak,  as  shown  in  Figure  7, 
from  ~2300  MLT  to  ~1 800  MLT. 

[55]  4.  Using  the  Siscoe  [1982]  ring  current  model  (S-RC) 
and  knowing  the  time-dependence  of  J1  and  J2,  a  time 
profile  for  the  penetration  electric  field  Eeq(t)  can  be  found 
as  shown  in  Figure  9.  From  Figure  9  it  is  seen  that  the  earlier 
bubbles  are  associated  with  values  of  Eeq  >  1 .25  mV/m, 
consistent  with  the  Fejer-Whalen  threshold  [Fejer  et  al., 
1999;  Whalen,  2001],  although  it  should  be  emphasized 
that  other  mechanisms  exist,  such  as  ionospheric  winds 
[. Richmond  et  al.,  2003],  as  alternate  interpretations. 

[56]  We  have  compared  and  combined  the  Nopper-Car- 
ovillano  [N-C],  the  Hill-Siscoe-Ober,[H-S-0],  Siscoe  Ring 
Current  (S-RC),  and  Weimer  models  in  a  consistent  way 
(see  Figure  6).  In  this  manner,  the  solar  wind  has  been 
electrically  coupled  to  the  ionosphere. 

Appendix  A 

[57]  In  this  appendix  we  have  two  goals.  The  first  goal  is 
to  illustrate  some  of  the  global  features  of  N-C.  The  second 
goal  is  to  demonstrate  how  ionospheric  electric  fields  and 
currents  can  be  scaled  with  the  transpolar  potential  using 
Table  1  and  (2).  In  this  way,  estimates  may  be  made,  for 
example,  of  electrojet  intensity  as  a  function  of  solar  wind 
parameters.  Figure  Al  shows  the  electric  field  vector  pattern 
at  high  latitudes,  as  obtained  from  N-C.  Recall  for  this  case 
that  the  input  to  N-C  is  a  uniform  region- 1  current,  forming 
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Figure  A2.  The  ionospheric  current  in  amperes  per  meter 
at  high  latitudes.  Note  that  the  dominance  of  the  Hall  term 
that  leads  to  electrojet-type  properties  in  the  auroral  region, 
a  westward  electrojet  on  the  dawnside  with  an  eastward 
electrojet  on  the  duskside.  The  MLT  axis  has  been  expanded 
to  accommodate  the  arrows  at  the  figure’s  edges. 


a  longitudinal-arc  of  120°  and  extending  3°  in  latitude  into 
the  ionosphere  on  the  dawnside,  out  of  the  ionosphere  on 
the  duskside  (see  Figures  3  and  4).  The  N-C  model  clearly 
show  in  Figure  A1  that,  with  the  ionospheric  conductance 
given  in  (17),  this  is  equivalent  to  placing  a  positive  charge 
at  0600  MLT  and  74°  in  magnetic  latitude  and  a  negative 
charge  of  equal  magnitude  at  1900  MLT,  also  at  74°. 
Therefore  given  the  ionospheric  conductance  profile  accord¬ 


ing  to  N-C  there  is  a  correlation  between  the  geometrical 
distribution  of  the  region- 1  currents  and  the  two-cell  polar 
potential  pattern.  Returning  to  Figure  Al,  the  distance 
between  the  two  charges  is  30°  or  3.3  x  106  m.  If  one 
now  takes  the  value  of  73  kV  for  <E>PC  as  given  in 
Table  1  then  the  average  electric  field  across  the  polar  cap 
is  22  mV/m,  consistent  with  the  arrow  lengths  shown.  If, 
for  example,  one  wants  to  know  the  ionospheric  electric 
field  when  <E>/>C  =  180  kV,  as  determined  from  (2),  one 
finds  from  Table  1,  J1  =  (13.7  x  0.180)  =  2.5  MA.  This 
changes  the  scale  factor  in  Figure  Al  from  20  mV/m  to 
50  mV/m  and  the  average  electric  field  strength  to  55  mV/m. 
In  this  way,  the  ionospheric  electric  fields  can  be  estimated 
as  a  function  of  the  solar  wind  parameters  to  a  zero-order 
approximation.  The  tenn  zero-order  refers  to  assuming  the 
polar  cap  conductance  being  static  throughout  the  storm 
and  the  neglect  of  other  nonlinear  processes. 

[58]  Now  let  us  look  at  the  high-latitude  currents.  The 
ionospheric  Hall  conductance,  when  impressed  by  the 
electric  field  shown  in  Figure  Al,  causes  the  ionospheric 
currents  to  form  the  pattern  shown  in  Figure  A2,  a  westward 
electrojet  on  the  dawnside  and  an  eastward  electrojet  on  the 
duskside.  The  currents  are  predominantly  Hall  currents, 
forming  closed  loops  over  the  poles  and  at  lower  latitudes. 
Under  the  assumption  that  the  electrojets  are  concentrated 
between  60°  and  72°  in  magnetic  latitude,  Figure  A2 
implies  an  eastward  electrojet  of  1.0  MA  compared  with 
0.6  MA  for  the  westward  electrojet.  On  the  other  hand,  if 
one  now  numerically  integrates  the  azimuthal  current  den¬ 
sity  from  0°  to  72°  degrees,  one  finds  a  total  eastward 
current  of  1 .4  MA  at  dusk  and  a  total  westward  current  of 
1.1  MA  at  dawn.  In  other  words,  in  this  model  71%  of  the 
eastward  current  exists  above  60°  in  magnetic  latitude  in 
comparison  with  55%  for  the  westward  current.  If  we  now 
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Figure  A3,  (a)  The  latitudinal  Ex  and  azimuthal  components  of  the  electric  field  for  the  two  cases 
shown  at  1800  MLT.  Note  that  the  superposition  of  the  two  cases  shows  the  suppression  (shielding)  of  J1 
by  J2.  (b)  The  same  quantities  at  higher  latitudes.  The  observed  divergence  in  Ex  shows  that,  for  the 
simple  current  configuration  chosen,  the  precipitation  current  manifests  itself  as  a  point  charge,  leading  to 
the  potential  contours  seen  in  Figures  3  and  4.  The  negative  values  of  Ex  near  the  pole  arise  from  the 
transpolar  potential. 
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multiply  these  results  by  2.5  for  the  6-7  April  2000 
magnetic  storm,  as  was  done  above  for  the  electric  field, 
we  may  estimate  the  electrojet  current  values  at  maximum 
<I> PC.  These  values  are  2.5  MA  and  1.5  MA  for  the  eastward 
and  westward  electrojets,  respectively. 

[59]  Figure  A3  shows  the  electric  field  profile  along  the 
1800  MLT  meridian  as  predicted  by  N-C.  Figure  A3a  shows 
two  cases  with  either  J 1  or  J2  =  1  MA  with  the  other  being 
zero.  The  symbols  represent  the  actual  output  of  the  N-C 
code.  The  symbol  Ex  denotes  the  poleward  electric  field 
component  along  the  magnetic  meridian  and  E^  the 
eastward  component  which,  under  symmetry,  is  equal  to 
Eeq.  Nopper  and  Carovillano  [1978],  by  requiring 
symmetry,  imposed  the  boundary  condition  Ex  =  -Eg  =  0 
at  the  equator  and  solved  for  one  hemisphere.  We,  on  the 
other  hand,  retain  symmetry  but  treat  both  hemispheres  and 
do  not  impose  an  equatorial  boundary  condition.  Note, 
however,  that  for  both  cases  in  Figure  A3a,  Ex  — >  0  at  the 
equator  (Xm  =  0)  consistent  with  symmetry. 

[60]  Now  one  of  the  main  goals  of  this  paper  is  to 
estimate  Eeq  under  various  high-latitude  conditions, 
where  Eeq  is  the  difference  between  the  equatorial  E9 
due  to  J1  and  J2.  If  Eeq  is  positive,  it  contributes  to 
stimulating  the  Rayleigh-Taylor  instability  for  bubble 
formation.  On  the  other  hand,  if  Eeq  is  negative  it 
suppresses  the  instability  (overshielding).  Figure  Ala  shows 
that  for  J1  =  J2  =  1  MA  that  Eeq  =  -  0.4  mV/m,  which  is 
consistent  with  overshielding. 

[61]  Figure  A3b  shows  the  high-latitude  electric  field 
profiles  for  J1  =  1  MA,  J2  =  0  MA.  Recall  that  for  the 
N-C  model  in  the  present  instance,  the  region- 1  currents  are 
assumed  to  be  uniform  and  form  a  semicircular  slab  about 
the  dawn  and  dusk  meridians.  Under  this  simplistic  condi¬ 
tion  we  see  from  Figure  A3b  that  there  is  an  equivalent 
charge  density  located  at  the  center  of  the  slab.  The  rather 
dramatic  discontinuity  in  the  slope  of  Ex  is  unrealistic  as  a 
result  of  this  simplistic  assumption  and  can  easily  be 
removed  by  adjusting  the  distribution  of  Jl.  Poleward  Ex 
becomes  rapidly  negative  to  reflect  the  dawn-dusk  potential 
across  the  polar  cap.  Note  that  Ex  dominates  over  E:p 
everywhere  but  at  the  ionospheric  equator  where,  as  just 
pointed  out,  Ex  goes  to  zero  because  of  symmetry. 

[62]  In  this  appendix  we  have  illustrated  that  N-C  despite 
its  conceptual  simplicity  is  a  powerful  tool  in  quickly 
defining  the  electrical  properties  of  the  global  ionosphere, 
as  modified  by  the  solar  wind.  It  is  argued  that  N-C  gives 
reasonable  zero-order  estimates  for  the  ionospheric  electric 
field  and  may  be  modified  in  a  systematic  way  to  explain 
detailed  observations. 
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